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Abstract: Index theorems for the Dirac operator allow one to study spinors on manifolds with 



is computed using a supersymmetric quantum mechanics representation. Here we find agreement 
with existing results which employed heat kernel and Pauli-Villars techniques. Nonetheless, this 
computation also provides a stringent check of the Feynman rules of de Boer et al. for the 
computation of quantum mechanical path integrals. Our results can be verified via a duality 
relation between manifolds admitting a Killing- Yano tensor and manifolds with torsion. As 
an explicit example, we compute the indices of Taub-NUT and its dual constructed using this 
method and find agreement for any finite radius to the boundary. We also suggest a resolution 
to the problematic appearance of the Nieh-Yan invariant multiplied by the regulator (mass)^ in 
computations of the chiral gravitational anomaly coupled to torsion. 
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1. Introduction 

Index theorems for manifolds with boundary and torsion involve a detailed understanding of 
spinors on such manifolds. In the torsion free case this problem was understood long ago in the 
mathematics literature and is solved by the Atiyah-Patodi-Singer index theorem Essentially, 
they found that the index for manifolds with boundary is the sum of the usual bulk result for 
the index (proportional to the integral over the first Pontrjagin class) plus a correction from 
the so-called i] invariant of the boundary Dirac operator. Upon the addition of torsion, further 
questions arise such as how does torsion affect the computation of the index in the bulk? What 
are the appropriate boundary conditions for spinors on such manifolds and how does one define 
the adjoint of operators acting on these spinors? How does torsion modify the computation of 
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the 7] invariant? Are there corrections to the index over and above that found by Gilkey if 
the metric is not of product form near the boundary? From a physical viewpoint the resolution 
of these questions is of considerable interest, not only because of the relation between the index 
of the Dirac operator and anomalies, but also due to the considerable importance attached to 
the physics of supergravity theories (involving, of course, both spinors and torsion) on manifolds 
with boundary, e.g. anti de Sitter spaces [Q. 

Manifolds in the Einstein-Cartan theory which include torsion are not nearly as well studied 
as those without, and solutions of the equations of motion tend to be of considerable complexity, 
making explicit checks of the general formalism we are about to present cumbersome. Fortu- 
nately, we can relate the index of a manifold A4 with torsion to a torsion free manifold A4 when 
their Dirac operators, p (e. A) and p (e, V) respectively, satisfy ^ 

{p{e,A),p{e,V)} = (1.1) 

(note that p (e, V) denotes the Dirac operator with inverse vierbein e.r'^ coupled minimally to 
an abelian vector field V^, whose necessity will be explained later, and p (e. A) refers to inverse 
vierbein and coupling to totally antisymmetric torsion which may be re-expressed via Hodge 
duality as minimal coupling to an axial vector field A^) . This equation will be made mathemati- 
cally more precise later, it suffices to say here that if we are able to find a pair of manifolds whose 
Dirac operators are related in this way, we will be able to test our results for index theorems on 
manifolds with boundary in a concrete example. In what follows we will refer to the manifold 
with torsion A4 as the "dual" of the torsion-less manifold A4. 

In fact, as was pointed out in P], for the case where one takes the manifold Ai as the 
background for classical spinning particle models with extended supersymmetry, one may view 
p (e, V) = Y^T^'^ii + ■ ■ ■ a supersymmetry generator and then search for additional super- 
symmetry generators p (e. A) = Y^r'^dfj^ + ■ ■ ■ . The tensor e^^^ is known as a Killing- Yano tensor 
and may be viewed as the inverse vierbein of the dual manifold Ai. We also stress that in the 
quantum case, when one considers the Dirac operators p (e. A) and p (e, V) depending on Dirac 
matrices satisfying a Clifford algebra, equation (|1.1|) is more subtle than in the classical case 
where one has supersymmetry charges depending on world line Majorana spinors obeying canon- 
ical Poisson brackets. The quantisation of Killing- Yano symmetries is, in fact, an important 
aspect of our work. 

At the end of this paper we take as an example to be Taub-NUT, 



ds 



_r, r + 2m 



dr^ + r^ de^ + r'^ sin^^d^^ 



Arm 

+ 



r + 2m 



dip + cos 9 d(f) 



(1.2) 



^The existence of a relation of this type was first suggested by Jan-Willcm van Holtcn based on the vanishing 
classical Poisson brackets between the supercharges of manifolds with Killing- Yano supersymmetries B . 
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a manifold of considerable physical interest and A4 is then the dual Taub-NUT manifold with 
metric 



r + 2m r , o m^r^ / , „o . o „ , , ox 1 4rm^ 



+ 7 (d^^ + sin2^d0^) + — d^ + cos^d0 , (1.3) 



r + 2m 



along with totally antisymmetric torsion (our normalisations can be found in the appendix) 

T = -^P^^deAd<f>Ad^^. (1.4) 
(r + m)^ 

We should mention that the dual Taub-NUT manifold possesses neither self-dual Riemann curva- 
ture tensors nor satisfies the Einstein-Cartan equations but this is unimportant for our purposes 
since it does satisfy the central equation ( |1.1|) . 

Our chief line of argument to compute the index of a manifold with torsion and boundary 
follows the Atiyah-Patodi-Singer index theorem |]l|. We split the computation into a contribution 
from the bulk and a boundary contribution. The bulk contribution is computed by representing 
the index as the regulated trace of 7^ 



index ^ = Tr 



7^exp(^^^) 



:i.5) 



which in turn may be represented as a supersymmetric quantum mechanical path integral ||^. 
In this approach, however, one essentially inserts infinitely many sets of plane wave resolutions 
of unity, for which one may clearly integrate by parts. Therefore this calculation is insensitive 
to the precise boundary conditions obeyed by the states appearing in the trace in ( |1.5| ) and 
yields only the bulk contribution to the index. In fact, analogous computations have also been 
carried out in a heat kernel as well as Pauli-Villars framework and yield results identical to 
ours. Nonetheless, our independent computation is interesting in its own right since it provides 
a stringent test of the rigorous quantum mechanical path integral formalism for sigma models 
in curved space developed by de Boer et al. [^, [T^ (and in particular of the Christ and Lee 
counter terms ||ll[ that must be added to the action in order that the path integral provides 



a faithful representation of the trace (p..5|)). The details of this computation are described in 
section |^. 

Having calculated the bulk contribution one must compute the boundary contribution which 
has two parts. The first is the rj invariant for the boundary manifold. The r] invariant is a 
spectral invariant of the boundary Dirac operator and may be computed from the eigenvalues of 
this operator. In the torsion free case this problem was solved for various asymptotic geometries 
by Hitchin (and applied to the Taub-NUT manifold in 0, |13|, Q)- Therefore, we generalise the 
computation of Hitchin to include the dual Taub-NUT example. To this end we construct and 
solve the boundary Dirac operator in the presence of torsion which allows us to compute the 77 
invariant for this manifold. 
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The second boundary contribution, due to Gilkey 0, is necessary in the case where the 
manifold is not of product form at the boundary. The main idea is that the APS analysis is valid 
for manifolds which approach a cylinder whose cross section is the boundary manifold. If this is 
not the case, one must make an extra boundary correction proportional to the integral over the 
boundary of the difference between the Chern-Simons form computed using the product metric 
and that computed using the metric itself. This computation can be easily generalised to the 
torsion-full case where we find a boundary correction proportional to a generalised Chern-Simons 
term. Note, in particular, that all our results hold for manifolds with boundaries at any finite 
radius, rather than just in some non-compact limit. The derivations of the boundary corrections 
are to be found in section 0. 

Section ^ analyses in detail the precise relation between index theorems of manifolds pos- 
sessing mutually anticommuting Dirac operators. This relation requires a careful understanding 
of boundary conditions for spinors and therefore the space of states in the trace ( p..5| ). Sub- 
section outlines the main ingredients of the Killing- Yano technique to generate manifolds 
with mutually anticommuting Dirac operators and the explicit geometric data for our example 
manifolds, Taub-NUT and dual Taub-NUT may be found in |4.3| . Also included in subsection O 
is the generalisation of the Killing- Yano construction to the quantum case. 

Urgent readers not wishing to wade through all the details may find the final expression for 
the index with torsion on manifolds with boundary in equation ( p.32| ). A summary of the results 
we obtained for the example manifolds is presented in subsection ^4.4 Finally, more speculative 
remarks along with possible applications of our results are reserved for the conclusion where, 
in particular, we discuss the controversial 1/(3 regulator (mass)^ Nieh-Yan contribution to the 
index in light of our work. In passing, we note that many of the issues we encounter have been 
dealt with in some depth in the mathematics and physics literature so that we have attempted 
to present these more formal aspects in an accessible format. 



2. Bulk contribution to the index with torsion 
2.1 The Dirac index as a quantum mechanical path integral 

This section is devoted to the computation of the bulk computation to the index employing the 
supersymmetric path integral approach of Alvarez-Gaume and Witten . The bulk contribution 
to the index of the Dirac operator with torsion has already been computed in via a heat 
kernel expansion and in using Pauli-Villars regularisation, needless to say, our results exactly 
reproduce those of these authors. Besides the high degree of certainty one now has in the result 
(see equation (|2.19|) ), our motivations for performing this calculation were fourfold: 

1. Our computation is carried out in the framework of rigorously defined quantum mechanical 
curved space path integrals introduced by de Boer et al. [1^ and therefore provides a 
stringent check of their work. In particular, we note that the result depends crucially on 
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the Christ and Lee terms [|Tl|] (see the second hne of (|2.13|) ) that one must add to the action 
appearing in the path integral at order = in order to employ propagators derived via 
a midpoint rule. 

There has been some controversy in the literature over the appearance of a term in the 
(bulk) index with torsion that diverges as one takes the limit in which the regulator is 



removed ||15| , |16|, [T^, [l^, |20|. This term is proportional to the Nieh-Yan tensor (see 
the first term on the right hand side of ( p.l9D ) and arises at tree level in our path integral 
expansion. In the path integral approach, one clearly sees the importance of including this 
term in order to obtain the correct results for the bulk index [i.e. the regulator independent 
contribution) since it appears as a vertex as well as in disconnected graphs. 

3. The path integral derivation is based on inserting infinitely many plane wave resolutions 
of unity. Therefore one clearly sees that it is only sensitive to the bulk contribution, since 
differentiations by parts then ignore any boundary terms. 

4. The quantum mechanical model itself is of intrinsic interest since it describes a spinning 
particle moving through a background with torsion. Obviously, the construction of the 
precise path integral quantisation of this model is an important step in understanding its 
quantum dynamics. 

The index of the Dirac operator p is given by the trace over the spectrum 



index ^ = Tr 



7^ exp (^-p '■ 



(2.1) 



When there is a gap between the zero-modes and the remainder of the spectrum |2T[] this result 
should, of course, be /3-independent. However, for manifolds with boundary, one must carefully 
take into account the states traced over. If one is interested in the bulk contribution only, the 
trace can be evaluated by inserting N sets of plane waves (where dt ~ — > 0) and the result 
may then be represented as a path integral. The expression one then finds is an expansion in 
powers of (3 

/?-'A_i(^) + Ao(^) + 0(/3) = A(^). (2.2) 

The one- loop beta independent contribution Ao{p) yields the bulk result for the index and as 
will be explained in section ^ the complete result for the index is obtained by including the APS 
and Chern-Simons boundary correction terms. Here we will compute the terms A_i{p) and 
Af){p). We comment on the interpretation of the 1/(3 contribution A_i(^) in the conclusion. 

Before introducing the details of our path integral computation, let us briefly summarise the 
main ingredients of the approach. To begin with, one treats the exponential in the trace ( |2.1| ) 
as the imaginary time evolution operator of a quantum mechanical system. The operator in 
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the exponent corresponds to the Hamiltonian. Importantly, its operator ordering is now fixed 
so that the usual step where one begins with some classical Hamiltonian generating dynamics 
via Poisson brackets and then elevates this expression to an operator, fixing the ordering by 
symmetry principles (or ultimately experiment!) is not needed. The operator H = —j3'^p^/2 is 
the quantum Hamiltonian. 

The next step is to derive the path integral by inserting N complete sets of states. The 
achievement of de Boer et al. was to derive a path integral which faithfully represents the operator 
expression The point being that once one takes the continuum limit — >■ oo, the discretised 

expressions derived at finite N become distributions whose products (which appear in Feynman 
diagrams) are ambiguous. However, given the rigorous discretised expression for the path integral 
it is always possible to resolve these ambiguities in a unique and finite fashion. Their derivation 
was performed by rewriting the quantum Hamiltonian in a Weyl ordered form which allows 
infinitesimal transition elements at intermediate steps in the path to be evaluated by a midpoint 
rule which leads to a certain set of vertices, propagators and rules for products of distributions. 
Of course another ordering principle would have yielded different propagators and vertices, the 
key point is to derive a path integral which precisely represents^ (|2.1|). 



We now turn to the details our calculation and begin by thinking of the exponential in (pTT 



expH(-/?2^V2)(-z)//?] (2.3) 

as a quantum mechanical evolution operator exp[— i/ft//?,] with Hamiltonian — /3^^^/2, unit 
Euclidean time interval and Planck's constant h = (3. The Dirac matrices are then identified 
with fermionic coordinates and the derivative operator with the canonical momentum of the 
bosonic coordinate x^, 

(2.4) 



^It is interesting to compare this situation to what one finds in quantum field theory. In quantum field 
theory, one must regulate and renormalise to obtain a path integral yielding finite, well defined results for the 
products of distributions appearing at loop level in Feynman diagrams. Furthermore, one makes some choice of 
renormalisation point in order to make contact with physical quantities. Of course, in quantum mechanics, there 



is no anomaly in scale invariance and the path integral Z[y^ ^ W] (see (2.6)) should provide an unambiguous and 
finite representation of the transition amplitude (y, ^'[7^ exp^^^ \y,'^). Although, the quantum mechanical 
Hamiltonian H = —0^Ip'^/2 has a definite operator ordering, in passing to a path integral representation, the 
vertices of the path integral are ambiguous up to reorderings of the Hamiltonian. However, once a particular 
ordering (there is nothing sacrosanct about Weyl ordering, it is merely convenient) is fixed one must decide 
(derive) which rule is used to define products of distributions. In analogy with quantum field theory theory, the 
discretisation procedure and rule for the products of delta and Heaviside functions above may be thought of as 
regulating and renormalising the theory with renormalisation point judiciously chosen to reproduce the required 
transition amplitude with the given Hamiltonian. 
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The similarity transformation for the canonical momentum (with g = det g^iy) implies the 
inner product {x\y) = g~^^'^5^{x — y). Observe also that the rescalings by appropriate factors of 
13 in ( p.4| ) ensure that (3 plays the role of Planck's constant in the canonical commutation relations 
and that the Hamiltonian —j3'^p'^/2 begins with a "classical" (3 independent contribution, plus 
possibly terms linear and quadratic in (3 (depending on the precise ordering in which one writes 
the operators x", and 

We may now represent the index by a path integral with periodic boundary conditions, 
the presence of 7^ in the trace yields periodic boundary conditions for the fermions also. 



A{p)= I dS d'^ Z[y", v&T (2.5) 

where the path integral Z is a function of constant real background fields y'^ and Schemati- 
cally 



Z[y>', ^^^] = j [dq>' dri' da'^ dfe^ dc^] 



expj-^ [ dt/:(x'^(t) =y^ + g'^(t),^"(t) = ^" + r/"(t))| 
^ P Jo ' 



(2.6) 



As discussed above, one usually expects path integrals to be plagued with ambiguities associated 
with the normalisation and precise definition of the measure along with products of distributions 
in loops. Fortunately however, finite interval path integrals in curved space have been studied in 
detail by de Boer et al. |P, |10|. In particular they have found the precise definition of Z above, 
i.e. the exact vertices and Feynman rules, such that the path integral representation (|2.5| ) is 
identically equal to the index A(^ ). A two loop verification of their work in the torsion-less case 
may be found in p2|] along with the conventions employed in this paper. We now spell out the 
key details of this approach. 

Firstly, although manifestly the quantum mechanical trace ( |2.1| ) is finite, in the path integral, 
closed bose loops {q"{t)(f{t)) yield delta function divergences 5(0) (first observed by Lee and 
Yang P3|). However, the decomposition of unity 1 = jd^x g^/'^\x) {x\ yields a measure factor 



g^/'^ at each point of the path which one may exponentiate via bosonic and fermionic ghosts a"{t) 
and {6^(t), c^(t)}, respectively ||2^, ^ yielding the ghost action 



1 .1 1 



^0 



dt/:ghost= / dt-g{x{t)),,{a{tra{ty + 2h{trc{tY) , (2.7) 



whose net effect is to precisely cancel the divergence 5(0) whenever it appears (for example graphs 
including a ghost loop are needed to cancel the divergence in the third graph of table H). 

Secondly, ambiguous products of distributions involving Heaviside and delta functions 5{t — 
s)6{t — s) also appear. The result of P, ^ is that adopting vertices corresponding to a Weyl 



ordered Hamiltonian, one finds propagators such that 6'(0) = 1/2 (Strictly speaking, one returns 
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to the discretised derivation of the path integral where the delta function is a Kronecker delta 
5ij and the propagators depend on the discrete Heaviside function with On = 1/2). More pre- 
cisely, one must rewrite the operator valued Hamiltonian H = —(3'^p'^/2 in Weyl ordered form. 
Thereafter one obtains the Lagrangian^ 

C = -ip^x>' + ^tlJrr + H. (2.8) 

by replacing operators by c-numbers in the Weyl ordered Hamiltonian. In general one then finds 
additional vertices, over and above the classical result, of order (3 and The exact propagators 
are then those given in |10[ and in our conventions in |2^. Note that employing the vertices 
of (|2.13| ) and propagators of table |l| along with the above rule for products of distributions even 



the overall normalisation of the path integral as given in ( p.5| ) is exactly correct. 

Our task then is to rewrite the operator H = ^/2 in Weyl ordered form, i.e. in terms 

of symmetrised products of the bose operators and and antisymmetrised products for the 
fermions ip^. In particular, for the fermions, this just means that all products of Dirac matrices 
should be written in the basis 7'', 7*"*, 7''** and 7^^*". This computation will involve three different 
spin connections so to avoid confusion let us spell out our notations. The Dirac operator is 

p =p{e,A)= YDiu), , 7^^ = e.V , Di^), = d, + loo.rsY' ■ (2.9) 

Here D{u)^ denotes an operator depending on the spin connection uj^^rs = i^(e)^rs + ^firs where 
'^(e)^rs is the torsion free connection and A^^rs is a completely antisymmetric contortion ten- 
sor. The classical (/3 independent) contribution to H has been computed in and yields the 
Hamiltonian for a spinning particle depending on a third connection cD^^s = ^^/irs + '^A^^-s = 
'^(e)^rs + 'iAfj^rs with contortion three times the usual one. Hence we introduce the operator 
D(Q)n identical to the previous definition except depending on the new connection u^rs- We 
reserve the symbol D^{= e^/D^) for the torsion free covariant derivative. 
We now make some simple algebraic manipulations 

p 72 = ]^{g^^^D{uj), - V\>^ + i^^D{uj), + ArJ^YllD^u), (2.10) 



-g^^{D{u),D{u),-T,pf'D{u) 



1 



+ - {R{u),,r-s - 2A%,Aprs)l^-Y' + -{D^A^^rs)!''' (2-11) 



with R^j_yrs the curvature built from the spin connection uj^rs (our conventions for the Riemann 
tensor in the presence of torsion are spelled out in the appendix). The remainder of the Weyl 
ordering computation for the first term in ( |2.11| ) follows the computation in |]10| and the remaining 



•^Note that the exponent of the path integral is then — (l//3)/£ so that if one continues the Euchdean time 
t ^ it back to real time the exponent becomes i(l//3) J{pij,x^ + {i / 2)il'r'4'^ — H), i.e. none other than the usual 
Legendre transformation relating Lagrangian and Hamiltonian. 
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terms require only the identity 'j^'^'jr 
the identifications in (UA 



J rs 



+ 47!^ 



Hence we find, making 



1 r 



Z Z "I 1 

Z z J W Z 



w 



o Z 



(2.12) 



where [. . denotes Weyl ordering. We may now take the Weyl ordered expression for H and 
replace operators by c-numbers. Integrating out the momenta one finds the Lagrangian 



'■stu 



(2.13) 



At this point the vertices of the theory are those defined by (|2.13|) and propagators along with 
products of distributions are those found according to the midpoint rule. Observe that the term 
at order /3 vanishes identically. For zero torsion, the above result has appeared in p^. Note 



also that the first two 0(/9^) terms in the second line of ( |2.13| ) are neither general coordinate 
(GC) nor local Lorentz (LL) invariant. From first principles however, we know that the result 
for the index must enjoy these invariances, since varying the background fields Cf/ and cj^rs with 
respect to GC and LL transformations, the operator p"^ varies into a commutator [p'^,G + L], 
{G and L are the generators of GC and LL transformations) which vanishes under the trace 
as [7^,6* + L] = = [7^,^^]. The resolution of this apparent dichotomy is simply that the 
propagators, based on the midpoint rule obtained by Weyl ordering, are also not GC or LL 
covariant. Nonetheless, as we shall see, the combination of these two "evils" will yield a GC and 
LL invariant result. 



2.2 Loop expansion of the path integral 

We are now ready to compute the path integral via a loopwise expansion in the Planck constant (3 
to orders f3~^ and f3^ employing the Feynman rules for the finite-time path integral of |1^. The 
split into free and interacting parts is made via the expansion around backgrounds and \E' so 
that x^ + and ip"^ \E'''' + rj"^. The spin connection Q is now an independent background 
field hence in our diagrammatic computation we make a Riemann normal coordinate expansion 
for the metric only (details can be found in the references just mentioned as well as in [^]). We 
have listed the propagators in table |l] and some of the required integrals in table ||. Other useful 
diagrammatic tricks and identities for this type of computation have been developed in p2| . 

We only need compute graphs quartic in background fermions (in the diagrams they 
are denoted by a dash on the external line) in order to saturate the Grassmann integration d^\E' 
in (|2.5|) which in turn yields a four-dimensional Levi-Civita symbol. There are several "tree level" 
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contributions^ (those contributions arising from the terms in the action that are independent of 
quantum fields). To order 



^tree 



M 



lb 



(2.14) 



beginning with the leading l/f3 contribution to which we draw the reader's attention. This 
contribution is proportional to the so-called Nieh-Yan tensor As expressed above it is 

clear that this tensor vanishes for exact torsion; alternatively one can re-express the Nieh-Yan 
contribution into a perhaps more familiar form by writing 



2di„A 



2 A\]^^yA^ pfj-^ . 



(2.15) 



where R{uj)pyp„ is the Riemann curvature of the torsion-full connection oj^rs = ^{.^)prs + A^^s- 
We will discuss this contribution, which has led to some confusion in the literature, in more 
detail in the conclusion. Observe that the order 13^ terms are also proportional to the curl of 
the torsion, though not immediately recognisable as a topological invariant as they are only a 
part of the result at this order in the Planck constant. Furthermore, observe that the product 
of the Nieh-Yan tree level tree level j3~^ contribution and the Christ and Lee Weyl ordering 
counter terms yield the second and third 0{j3^) terms in ( p.l4| ). This provides the long sought 
after explicit confirmation that these two loop counter terms (they appear with a factor in 
the action ( |2.13| )) are actually necessary in a rigorous definition of the path integral. 

The remaining graphs involve one- and two-loop integrals and are given in tables ^ and ^ 
Note that we have not drawn graphs obtained by different locations of the dots, which denote 
time derivatives on propagators, although these contributions are, of course, included in the 
quoted results. Adding all these contributions together our result for the regulated trace of 75 is 



M 



V9 
(27rz)' 



^rstu 



1 1 

J^DrAgt^ + — 



pv 



tu 



+ 8 Df'DpDrA 



■stu 



+ 36 A 



A""""^ DrA 



stu 



+ ARDrA 



stu 



+ 0{p) . (2.16) 



Note that the contracted indices in the term quadratic in the Riemann tensors are curved ones 
since in the presence of torsion the R{u})p,i,rs 7^ R{yj)rsp.v (see ( |A.15|) ). The graphs for this 



term are not displayed in the table but can be easily computed along the lines of [^]. Also, 

"^To avoid confusion, note that we only compute contributions of order and Z?", however, since we must 
compute the partition function Z[y^^ ^f^], the generating function of all diagrams, disconnected products of "two 
loop" and "tree" graphs will produce [3 independent results which must not be neglected. 
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observe that the scalar curvature on the second hne of (|2.16| ) is that of the torsion-free Riemann 
curvature since it comes from the figure eight graph on the last line of table ^ in which one 
expands the metric gi_iu{y + q) to second order in quantum fluctuations q{t)^ which in Riemann 
normal coordinates yields the torsion free Riemann curvature. 

A useful test of the result ( |2.16| ) is independence of the choice of complexification of the 
original Majorana spinors required to employ a coherent state formalism for the fermions. The 
three terms from the one-loop diagrams and the single term from the two-loop diagram do indeed 
vanish upon using the Schouten identity, by virtue of the fact that K"^"' is anti-symmetric. 

To facilitate comparison with existing results and see that the result is a total derivative it is 
useful to extract the dependence on the torsion free curvature which we denote by R^up'^ where 
Rpup^Va = [D^, Dy\Vp. Furthermore, we rewrite the result in terms of the axial contortion vector 



= ^rstu^'^*^- Using ( |A.15| ) one derives that 



R{^) pursR{^) tu fc 



rstu 



p T3 pa 

J^pupaJ^KX 



2D.A,D^Ay 



D.A.A'^A'' 



^PukX 



(2.17) 



The cross terms between R and the torsion vanish identically by virtue of the Ricci symmetry 
Rplupa] = and so do the terms quartic in the torsion. The second hne of (|2.16| ) is easily written 
as a covariant derivative by commuting the triple derivatives using 

D^'D^D^A" = D^D^'D^A'' + D^A^R^"" + lA^'D^R . (2.18) 

This brings ( p. 161) into a form derived previously by Obukhov using heat kernel methods^. 



M 



(27rz)2 {2/3 



—.D^A^^ 



1 1 
192 



^pCfKX 



RpupaR^ kX 



-F F \ 

pa^ kX 



1 

24 



+ — D^IC^' 



with 



[d'D. + IA'A. + IR^A^' , F^, = 2D^^A, 



+ 0{I3). (2.19) 



(2.20) 



In the case that the Nieh-Yan tensor vanishes, which is manifestly so when the torsion is 
exact, the above expression (|2.16| ) reduces to the term quadratic in R{uj) which coincides with 



the result derived by Mavromatos some time ago |26 



^Obukhov's work has been re-examined by Yajima and Kimura g^, ^ who found, again in a heat kernel 
approach, the same result for the bulk index with torsion but disagree upon the precise details of his derivation of 
the heat kernel coefficients and claim that it leads to incorrect results for the Lorentz anomaly in the presence of 
torsion. We do not wish to enter this debate here and it suffices to say that the path integral techniques presented 
in this paper can equally well be employed for other anomalies and could, therefore, provide an independent 
resolution of this issue. 
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M 1/ 



( qisTqity ) = g^-" t{l - s)9{s -t) + s{l - t)9{t - s) 



{ 7]{s)'^7]{tY ) = It]-^"" e{s -t)- 9{t - S) 



Table 1: Propagators for the action ( 2.13| ). The symbol K^' 
complexification of the Majorana spinors, as explained in 



parametrises the choice of 



As long as the K,^ vector is smooth (as it is for smooth metric and torsion), the additional 
contribution D^JC^ to the bulk index in the presence of torsion will only receive contributions 
from the boundary. The subtle interplay between the three order terms in (|2.19| ) can be 



clearly seen in the explicit example given in subsection 4.4.1 below 



We stress that the Nieh-Yan term in the action ( 2.13 ) along with the Christ and Lee two- 
loop Weyl ordering counter terms have been crucial in obtaining this result (see the two-loop 
diagrams in table ^ as well as ( |2.14| )). Although the 1/(3 Nieh-Yan bulk contribution to index 
has appeared in the literature before |jl5|, |l^, 0, |l^, |19|, in our path integral approach one 
sees clearly its importance for the order (3^ terms and hence the necessity to include it when 
computing the bulk index in the presence of torsion. We will comment further on this term in 
the conclusion. 

Finally, we observe that the folkloric statement; "the chiral anomaly is topological and 
therefore insensitive to the subtleties of the precise definition of the path integral" can be made 
very definitive in the approach employed above. The reason being that in the absence of torsion 
the different counter terms corresponding to varying ordering schemes have no affect on the 
result since they are all higher loop effects. Only once the torsion is turned on, so that 
the Lagrangian includes the 1//3 Nieh-Yan contribution must one worry about the subtleties of 
ordering schemes etc. when computing the (3 independent anomaly. 



3. Boundary contributions to the index with torsion 

3.1 The APS boundary correction with torsion 

On manifolds with boundary, the central question is which boundary conditions one must impose 
for spinors in order that the index problem is well posed. Exactly this problem was solved by 
Atiyah, Patodi and Singer in the mid-seventies ^ for manifolds which have a product form 
near the boundary {i.e. when the manifold is a cylinder with the boundary manifold as cross 
section). In this case the Hilbert space is sufficiently under control and consistent boundary 
conditions can be imposed on spinors. Moreover, due to the fact that the radial factor can 
be written down explicitly, the index computation can be reduced to the problem of suitably 
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n P Q 
n P Q 




Table 2: A number of useful integrals obtained by straightforward application of the prop- 
agators in table |l] where the definition of ^fC*"" may also be found. The remaining integrals 
required for our computation as well as a detailed explanation of the diagrammatic notation 
above can be found in |22]. 



counting the difference between the number of positive and negative eigenvalues of the boundary 
Dirac operator. 

For manifolds that are not of a product form, a split of the Hilbert space is not as obvious. 
Nevertheless, due to the topological nature of the index, it is possible to deform these manifolds 
until their boundary takes a product form. The index then splits into a bulk part plus the 
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n 1 ^ ^ rsnm-'- tu ~ -i rsnm-'- tupn-'^^ ; 

d4 Id 




Table 3: One-loop diagrams contributing to the path integral at order P^. Only those graphs 
that are not already zero because of symmetry properties of the background fields are displayed. 
Note that we denote = y^f'-^^y^f'^^''. 

APS boundary term, amended by a modification due to the deformation. This modification was 
pointed out by Gilkey [0] and we will generalise it here to also include torsion. 

Let us now explain these matters in more detail. Suppose d = d/dr is a vector normal to 
the boundary dAi of the manifold in question. Then the Dirac operator (possibly including a 
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(3/4) 
288 



rstu ) 





1 

192 



^ mnpq^^ ^ ± rstu 





4:RFrstu^ 



rstu 



Table 4: Disconnected two- loop diagrams contributing to the path integral at order (3^ . Note 

that the result for the diagram on the third line includes graphs obtained by permutation of the 

time derivatives. One must also add to this diagram the corresponding ghost diagrams which 

d2 

A(s,t) = (5(0) - 1. 



cancel the divergence 2X(s,s) 



dsdt 



coupling to the antisymmetric part of the contortion, assumed to have only components parallel 
to the boundary) takes the form 



|)(e,A)=7f^' (a+|a;(e)V]) +^ 



(3.1) 



where B evaluated at the boundary r = r;, is the boundary Dirac operator. We denote the curved 
index in the normal direction spanned by the coordinate r as [r] and uj{ey p[r] is the trace of the 
torsion-free part of the spin connection computed from the vierbein e^'^ (no confusion should arise 
between the flat index r and the radial coordinate r). Observe that since p (e, A) is antihermitean 
with respect to the volume element \/ g{e) the coefficient of the single outward pointing Dirac 
matrix 7M = Y^r^^^ is fixed to be that appearing in ( [3.11 ). However, to verify that the Dirac 
operator is antihermitean with respect to the inner product (S, \E') = Jj^df^x^^/gie) along 
with taking into account differentiations by parts that hit the measure factor \/ g{e) one must 
also ensure that any surface terms also vanish by imposing appropriate boundary conditions on 
spinors. However, local Neumann or Dirichlet boundary conditions do not lead to solutions of 
the first order homogeneous zero-mode Dirac equation. Instead non-local boundary conditions 
are required. 
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Let us write the Dirac operator ( |3.1|) in a chiral basis for the Dirac matrices 

where V is some function of the normal variable r. The 2x2 matrix B is hermitean acting on 
spinors on the compact boundary manifold dAi. The APS split of the index computation into 
a bulk piece and a boundary piece can only be achieved if (|3.2| ) takes a simpler form near the 
boundary. First of all, B has to be independent of r and the connection trace uj{e)'^p[r] should 
vanish. This holds when the manifold is a product metric near the boundary. In addition we will 
require that there is no "radial" component of the torsion [i.e. A[r]up = 0) as this would enter 
the Dirac operator in a way similar to the connection trace and prohibit an explicit solution of 
the radial part of the eigenfunction. Although the APS condition can be imposed by smoothly 
deforming the manifold near the boundary, it is not obvious that such a deformation argument 
can be used to set a radial torsion component to zero. 

When these conditions are satisfied, Atiyah et al. found that one can consistently impose the 
following non-local boundary conditions on spinors \E', 




> ^ = , (3.3) 



E 10(^1 0(r.)=O^ 

l>0 

E 10(^1 X(r.)=0 
z<o 

where the boundary is at r = and the two component spinors |/ ) are eigenspinors of B, 

B\l) = l\l). (3.4) 

Only for such non-local boundary conditions does a non-trivial set of zero-modes of the Dirac 
operator survive. The index splits in a bulk and a boundary piece as indicated according to 

index^ = index(bulk) - ^(^/^(O) + hs^ (3.5) 

where index(bulk) is the usual integral over the second Pontrjagin class, h is the number of 
harmonic spinors on the boundary manifold DM., i.e. the number of solutions to i? \l/ = 0, and 
?7b(s) is a spectral invariant of the boundary operator B given by 

sign(/) 

where the sum runs over the eigenvalues I of the boundary Dirac operator B. As a function of 
the complex variable s, ?7b(s) has a well defined analytic continuation to s = ^ 
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When B only has a finite number of eigenvalues one can take s 



before summing ( p.6|) (for 
). In this case, one 



other cases in which the APS conditions can be made explicit see e.g. p^ , ^ 
can understand the appearance of i]b{0) in the expression for the index rather easily. Assuming 
that near the boundary at r = rj, the manifold is a product metric, the zero mode equation on 
this cylinder can be written as 




Xk-l) Cfce'^nO. 



(3.7) 



(for the product metric we may simultaneously diagonalise B and d). This implies that positive 
chirality zero-modes (those for which only the upper components are non-vanishing) have k = I 
while those of negative chirality have k = The boundary conditions ( |3.3| ) imply that the 
positive chirality states have only negative eigenvalues of B while the negative chirality states 
carry positive -B-eigenvalues (due to the relation between k and / this implies that only zero- 
modes which vanish exponentially as r ^ oo remain). The contribution to the index from the 
cylinder is therefore the disparity between the number of negative and positive eigenvalues of B 
which is clearly given by (p.6|) at s = and this limit may be taken before performing the finite 
summation (the factor 1/2 in ( p.5|) arises from the fact that half of the spinor boundary modes 
do not continue smoothly into the interior 0, 0). The APS theorem renders this argument 
rigorous for operators B with an infinite number of eigenvalues and also shows that the index 
problem indeed splits in a bulk and boundary piece according to ( |3.5| ) by appropriately gluing 
the two pieces of the manifold. 

Before we discuss the Gilkey correction term, we will now first show how the rj invariant can 
be computed for a class of boundary Dirac operators general enough to include the torsion-full 
examples to be discussed later. 

3.2 Computation of the 7] invariant for manifolds with torsion and boundary 

The 7] invariant must be computed on a case by case basis and since it involves all non-zero 
boundary eigenvalues this can be rather tricky. Fortunately this computation has been carried 
out by Hitchin for a large class of boundary geometries |Q (to be specific, squashed metrics) . 
The addition of torsion modifies the boundary Dirac operator by what amounts to only a slight 
generalisation of Hitchin's calculations. However, in keeping with the rest of this paper and 
because the explicit solutions we find to the boundary Dirac equation will be needed when we 
explain how one is able to compare index problems for manifolds with mutually anticommuting 
Dirac operators in section ^, we will now present these additional details. 
We begin with a boundary operator of the form 



B = S (7=^^X1 + 7''S2) 



-7 



^3 



7 



1234 



4A 



+ Z 



(3. 
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where A, S and Z are arbitrary parameters (depending, in general, on rb)- The operators Sj 
commute with the S0{3) Killing vectors {i.e. they are the inverses of the left-invariant forms) 
and satisfy the algebra 

[Si,S2] = S3 and cyclic. (3.9) 

To motivate the somewhat obscure looking form we have chosen for the operator ( p.8|) , we note 
that it is exactly the case solved by Hitchin |3^ when S = Z = 0. The generalised form presented 
here can be handled too and (as we will see) appears in the example manifolds discussed in 
section ^. A useful basis for the Dirac matrices is 

y=f° M 7^=r°-^"^i f=r°-^"^i 7^=r°-^"0 (3.10 



10/ V icr.T / \ icTy / \ icTz 



SH^'y „°.,J. (3.11) 



One then has 

-P{-S,Z) 
where the "Hitchin-operator" is given by 

PiS, Z) = z^pia, + Esa,] + ^Sga, + + Z . (3.12) 

Therefore it is most useful to solve the eigenvalue problem P{S, Z) (f = zip. To begin with call 
Sj = —iJi where Jj is the angular momentum operator {i = 1, ..,3) (and we employ the usual 
eigenstates \j,m) of and J3). Note that the Hitchin operator commutes with the Casimir 
[J"^, P{S, Z)] = so we may look for solutions at a fixed j. Actually, it is easy to solve this 
equation with the following ansatz 

V^=f (3.13) 



I j, m+1) 

where we define \j,j + 1) = = |j, — j — 1). Then if m = j or m = — j — 1 we clearly obtain 
solutions with eigenvalue 

z = j/\ + (A' + 2)/4A + Z. (3.14) 



Now if j > m > — j — 1 so that A = (j — m)(j + m + 1) is real and non- vanishing then we 
have solutions with eigenvalues 

z± = a±SA~ {m + 1)/A + (A^ + 2)/4A + Z (3.15) 

where a± solves the equation 

2772 + 1 

a±-^;^«±-l = 0, (3.16) 
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so that explicitly 



z± — + Z . (3.17) 



To conform with the notation of Hitchin we relabel j = {p + q — l)/2 and m = (p — g — l)/2 so 
that the eigenvectors, eigenvalues and their multiplicities read 

eigenvector eigenvalue multiplicity 

1^'^)^ f ' ^ A + ^ + z 



I 2 ' 2 / / 



(3.18) 



where the integers p and g take all values strictly greater than zero (note that the multiplicities 
above refer to distinct eigenvalues). To understand the multiplicities quoted in ( p.l8| ), it is useful 
to briefly explain the representation theory of states on S^. The metric on the squashed 
(considered in this section) may be represented (up to an overall irrelevant scaling) in terms of 
left invariant forms {i = 1, .., 3) as 

ds'\^, = (Si)2 + (S2)2 + (SV(SA))2 (3.19) 

From the invariant forms one may write down left invariant vectors Sj which commute with 
the Killing vectors ki and satisfy the algebra (the explicit forms of these vectors are given in 
section |47^ ) 

[Ei,E2] = S3 and cyclic, 

[S,,A;,]=0, (3.20) 
[ki , ^2] = ks and cyclic. 

Representing the operators Sj on the space { \j, m) } as above and the Killing vectors on a second 
copy of this space { | J, M) } then from the relation 

Y,m' = j^ihf (3.21) 

i=l i=l 

expressing equality of the Casimirs on (which may be easily verified from the explicit forms 
in [1.3| ), one finds the representation space of functions on to be {|j, m) ® |j, M)} which yields 
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the multiplicities quoted above (since in ( p.l8| ) we suppressed the dependence on the "right 
representation" | J, M) ) . 

We are now left with the task of computing ?7_b(0) = ?7p(0) as in ( p.6| ) given the spec- 
trum ( p.l8| ). Therefore we must find the analytic continuation to s = of (we drop an overall 
factor (2A)'' which clearly does not modify ?7_b(0)) 

Vb{s) = — 5- 

^(p + AV2 + 2ZA)^ 



p + q 



p + q 



P,g>o 1^ (^A2/2 + ^D{XS,p,q) + 2Za) " (-AV2 + y/D{XS,p,q) - 2Za) ' 

(3.22) 

where D{XS,p,q) = ipq^XSY + {p — qY- (Note that in all cases that we are interested, A, 5* 
and Z are such that —X'^/2 + D{XS,p, q) — 2ZX > 0, which motivates the above split in a 
sum over positive and negative eigenvalues, but in general one should of course worry about 
this separation as well as any zero-modes of the boundary operator). The single sum can be 
expressed in zeta functions without any difficulty. The double sum requires a lengthy analysis, 
most easily performed by first Taylor expanding in powers of the new variable A defined by 
AV2 = AV2 + 2ZX, 

2s(s + l)(s + 2) /A2^^ 



2^ (y)/((^ + 1)/2) 



3! 



-j /((. + 3)/2) 



(3.23) 



where the dots vanish for s = and the function f{s) is defined by 

p + q 



p,q>0 



DiXS,p,qr- 



(3.24) 



The terms involving the zeta functions yield A^/4 — 1/6 while the analysis of the residues of f{s) 
is rather nontrivial and produces 

1 



res 



res 



/ 



/ 



s + 1 
2 

s + 3 



{{xsr - 1) , 



This finally yields the r] invariant at s = 



VBiO) 



1 

6 ~ 

A3 
3 



A 

3 

3Z 



16 Z'^ 



s=0 



s=0 



AZ 



(3.25) 



(A^)2 ■ 



3 



+ AS'Z-6Z 



12 Z^ 



A^ 

T 



12 Z' 
AS' 



+ 1 



352 



(3.26) 
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At 5 = 1 and Z = 0, we recover Hitchin's result //^(O) = —1/6 + — So summarising, 

our final result for the r] invariant for the class of Dirac operators in (|3.8|) is given by the above 



result (ICTp 



3.3 Boundary correction for non-product metrics 

The APS analysis presented in the previous two subsections is valid for the manifolds whose 
metric takes a product form 0] at the boundary (more intuitively, those manifolds which, near 
the boundary, are a cylinder with the boundary manifold as cross section). If this is not the case, 
the manifold has to be deformed to such a product structure near the boundary before the above 
machinery can be applied. The bulk term should in turn also be computed for this deformed 
metric. As an explicit form of the deformed metric is often not available, it is, however, much 
easier to determine instead the error that one has made by computing the bulk contribution 
using the original metric. This error term, due to Gilkey [0 for torsion-less manifolds, has to be 
subtracted from the bulk. 

The form of this correction is rather simple. In the absence of torsion, the bulk term can be 
written as a boundary integral 

— ^ [ RAR = —^ [ dC = —^ I C, (3.27) 

provided discontinuities of C are taken into account properly. For the metric deformed to a 
product structure near the boundary, we get a boundary integral with C replaced by the Chern- 
Simons term of the product metric f/^^, = g^y{r = rh), CK Therefore, in the absence of torsion, 
the correction term that has to be added is 

1 

24.87r2 

where the Chern-Simons three-form C = oj /\ R ~ ^uj /\ uj /\ uj (and similarly for C" built from 



[C» - C] , (3.28) 

dM 



"sharped" objects computed from the product metric). Note that the integrand of (|3.28|) can 



be written as 9 A R where 9 = u — u'^ is the second fundamental form and u and are the 
torsion-free spin connections computed from the metric and product metrics respectively. (Where 
we warn the reader that from here on we replace uj{e) by simply uj to denote the torsion- free 
spin connection.) However, for intuitive, along with practical reasons, we prefer the form given 

iH). 



m 



For the torsion-full case the integrand of the bulk index derived in section ^ above may be 
written as the exterior derivative of a generalised Chern-Simons form 

ujAR-^uAujAuj-^AAF-2}C = C{A) - 2 /C (3.29) 

where the three- form /C is the Hodge dual of the vector IC^ in (|2.20|) , 

)C>' = ^e^'''''''}C,p.. (3.30) 
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Therefore, the boundary correction for the torsion-full case is simply 

/ [C\A) - 2 /C» - C{A) + 2 /C] , (3.31) 
where K,^ is computed by inserting r = in the axial contortion vector before computing 



the covariant derivatives in (|2.2CI|) . 

3.4 Generalised APS index theorem 



Orchestrating the above bulk results presented in (p.l9|) as well as the t] term and the boundary 
contribution (|3.31|) we obtain an APS index theorem generalised to manifolds with torsion 



index^ (e, A) 



24.87r2 



M 



-F{A)AF{A)-2y^D,}C' 



+ 



24.87r2 



dM 



C«(A)-2/C»-C(A) + 2/C 



(3.32) 



with the three-form /C defined by (|3.3CI|) and (|2.2CI|) above. A result similar to ( |3.32|) has appeared 
in a string theoretical context in |^ for which one encounters exact torsion and most of the 
terms above vanish. Needless to say, in the general case presented here there are many extra 
subtleties as discussed in the preceeding text. 



4. Index theorems for manifolds with mutually anticommuting Dirac 
operators 

4.1 General formulation 

The existence of mutually anticommuting Dirac operators on a pair of manifolds can be employed 
to derive a relation between the index theorems on these manifolds which we shall now present in 
detail. At the end of this section we test this relation on the explicit example of the Taub-NUT 
and dual Taub-NUT manifolds. 

To begin with, given a manifold Ai with (possibly torsion-full) Dirac operator p (e. A) 
satisfying 

p (e. A) = -p (e. Ay , {p (e, A),^'} = 0, (4.1) 

then let us assume that we have found a second operator p (e, V) acting on spinors defined on 
JH which solves the equations 

{p{e,A),p{e,V)} = (4.2) 
p (e, V) = -p (e, Vy , {p (e, V), 7^} = . (4.3) 



22 



At this point we do not yet identify p{e,V) with the Dirac operator on a manifold Ai. Im- 
portantly the hermiticity requirement in ([4 .31) is formulated on the manifold A4 (the precise 
definition of the adjoint operation for spinors on manifolds with boundary was presented in sec- 
tion |3.1| ). In general, given that the inverse vierbeine e^^ 7^ Cr^, observe that it will certainly be 
necessary to include the additional vector coupling in p (e, V) to satisfy both ( |4.2| ) and ([4.3|) . 

As the operators —ip (e, V) and —ip (e, A) anticommute, we can construct a new operator 
'j^p (e, A) commuting with the original Dirac operator, 

[-tpie,V),j'Pie,A)]=0. (4.4) 

These two operators (hermitean on Ai) can thus be diagonalised simultaneously. Eigenspinors 
with non-zero eigenvalues occur in pairs for both operators, since they anticommute with 7^. 
However, eigenspinors with vanishing —ip (e, V) eigenvalues are not necessarily those with zero 
j^p (e. A) eigenvalues. Nonetheless, in the computation of the index, 

index p{e,V) = n+{p {B,V)) - n^{p {e,V)) , (4.5) 

they still occur in pairs (?t,+ denotes the number of zero- modes with positive chirality). The 
above expression therefore only receives contributions from those zero-modes that are also zero- 
modes of 'j^p (e. A), that is Ker' {—ip (e, V^)) C Ker' {•y^p (e. A)), where Ker' denotes the set 
of zero-modes which do not have a partner of opposite chirality. The same holds true when we 
interchange the role of the two Dirac operators and we therefore also have Ker' {—ip (e, V)) D 
Ker' {■y^p (e. A)). We may therefore conclude that the indices of —ip (e, V) and y^p (e. A) are 
equal and hence^ 

index p (g, V) = index p (e. A) , (4.6) 

since if \1/ is an eigenspinor of the operator 7^ p (e. A) then 7^\1/ is the corresponding eigenspinor 
oi-ip{e,A). 

The equation ([4.6[) becomes much more interesting if one is able to view the operator p (e, V) 



as the Dirac operator on the manifold Ai but one must then carefully consider the space on which 
the operators p (e. A) and p (e, V) act. The reason being that typically one is interesting in 
computing the index of the Dirac operator as a trace over states 



index ^ = Tr 



7' exp (^^p 



(4.7) 



where the trace extends over the space of spinors living on the manifold of interest. However, we 
do not want to compute the index of the second Dirac operator p (e, V) on the manifold Ai, but 
rather on the new manifold Ai obtained by expressing p (e, V) = Y Cr'^ 9^ + • • • and regarding 

^The above argument was pointed out to us by Jan-Willem van Holten M. 
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as the inverse vierbein. Only in this way do we obtain a relation between index theorems 
on independent manifolds A4 and Ai. A priori, however, there is no reason to expect that the 
space of spinors on Ai should coincide with that on Ai. Furthermore, although the operator 
p (e, V) was assumed to be antihermitean on the manifold Ai, since the volume elements of the 
two manifolds will in general not be equal, the operator p (e, V) will have no definite hermiticity 
viewed as an operator on spinors on At. Therefore we make the additional assumption, (which 
holds for the example we have in mind) 



Assumption: The space of spinors (which were defined in detail in subsection |3.1| in such a way 
that index problems are well-posed) viewed as a four-component space of functions on At and 
Ai coincide. 

Notice that as a consequence the index problem for the operator p (e, V) is well posed on the 
manifold Ai. This can be seen as follows. Even though the operator p{e,V) has no definite 
hermiticity on the manifold A4 the operator p (e, V) = \J g{e) / g{e) p (e, V) \J g{e) / g{e) is anti- 
hermitean on Ai. Furthermore if are a complete orthonormal set of eigenfunctions of p (e, V) 
on Ai then g{e) / g{e)'^ \ are a complete set of orthonormal eigenfunctions of p (e, V) on Ai. 
Clearly however, by the assumption, if the above similarity transformation is non-singular then 
the set of functions are still complete (although no longer orthonormal) on Ai. so that p (e, V) 
can still be diagonalised on Ai leading to a well posed index problem on that manifold. 

It should now be clear that given the above assumption we have equality of well posed index 
problems calculated on Ai and Ai, 

index ^ (e. A) = Tr^ 7^ exp {^p (e, Af) = Tr^ 7^ exp {^p (g, Vf) = index |) (g, V) (4.8) 

which is the main result of this section. If we recall the discussion of boundary conditions for 
spinors on manifolds with boundary in section |3.1| , we see that to verify our assumption for the 
manifolds Ai and Ai one only needs to require equality of the projection operators 

El^)(^l = (4.9) 

where |/) and |/ ) are the eigenvalues of the respective boundary Dirac operators B and B. 
Remarkably, we find that (|4.9|) holds for the Taub-NUT and dual Taub-NUT manifolds. 

Note that to compute the index of p (e, V) on Ai we can split the operator into a sum of 
antihermitean and hermitean operators on Ai by writing 

p{e,V) = p{e,0) (4.10) 

with 

p (e, 0) = [p (e, V) - \y{uj%, - u;{eY (4-11) 

V, = \^%,-\^{e)%, (4.12) 
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where uj{e)f^rs is the torsion free part of the spin connection computed from the inverse vierbein 
Cr'^. The trick now is that the hermitean term 7^^^^ is nothing but the couphng of the Dirac 
operator to a purely imaginary abehan gauge field whose contribution to the index we may obtain 
by a naive analytic continuation from the well known results for the index of the Dirac operator 
with such couplings^. 

In the remainder of this section we show that manifolds with mutually anticommuting Dirac 
operators can be found if invertible Killing- Yano tensor exist and the above relation is realised 
explicitly by Taub-NUT and its dual manifold. Also provided are all relevant geometric data for 
these manifolds. 

4.2 Killing— Yano dual manifolds 

Manifolds with torsion admitting two "Dirac operators" that satisfy the properties discussed 
in the previous section do in fact exist and are understood in a systematic way (at least 
in the case where the Dirac operator is viewed as a classical supercharge). Before we move on 
to an explicit index computation we describe the Killing- Yano technique for generating these 
manifolds. 

The operator p (e, V) can be viewed as the quantum analogue of the supercharge of a 
spinning particle and in this context the additional Dirac operator p (e, A) generates the extended 
supersymmetry. Up to ordering ambiguities this means that the conditions on the tensors 6// 
and Arst appearing in the additional Dirac operator such that it anticommutes with the original 
Dirac operator can be deduced from the existence conditions of extended supersymmetry of the 
classical model pil[j^(see for the extension to include torsion). In short, additional supercharges 
on the manifold Ai are given by 

Pie, A) = 7'-e/(a^ + \{oo{e),st - |(e-^)/c,,,)7^*) , (4-13) 
where the tensor Cr'^ satisfies 

D(e)^e/ + D(e),e/ = 0, (4.14) 

equation for e^^ on the manifold Ai. Here we have raised and 
lowered indices with the vierbein e^*^ on Ai so that e^,''" = e^j^Cs^e^'' and -D(e)^ denotes the 
torsion- free covariant derivative on M.. The connection is determined by requiring 

Crst = D{e)[rest\, (4.15) 

^In the example below we will anyway find a vanishing contribution from this additional abelian vector coupling, 
however in general it is quite simple to handle. The bulk term is just the usual result for the abelian anomaly 
modulo an overall minus sign due to the analytic continuation. The Chern-Simons and APS boundary corrections 
are computed exactly as described in section 0. 



which is the Killing- Yano |40 
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where again we have employed the vierbeine on Ai to flatten indices. Note however that when 
we come to regard Cr'^ as the inverse vierbein on the dual manifold A4 we will denote its inverse 
{i.e. the vierbein itself on A4) as e^^ = (e^"'^)^'". Viewing uj{e)^st — l{€.~^)f/Cpst as the connection 
on the dual manifold Ai the contortion on that manifold is then 

Kt^s = t^(e)^st - ^^{e)^,t - ^{e~^)/cpst ■ (4.16) 

the totally antisymmetric part (on Ai) of which defines the tensor Arst = K^rst]- 

Indeed just as complex structures (with one fiat index) can be used to define a local Lorentz 
frame different from the frame spanned by the inverse vierbein, an interpretation of the Killing- 
Yano tensor as an inverse vierbein is possible as well. In this case, however, the manifold for 
which the inverse vierbein is Cf/ is not identical to the original one given by e^'^ (the square e^^^e^r 
is a Killing tensor which does not normally coincide with the metric on A^). This observation 
was first elucidated in and plays a central role in this paper by providing an example of 
manifolds with mutually anticommuting Dirac operators. 

At the quantum level one must study the possible orderings when making the transition from 
a classical supercharge to the Dirac operator acting on a spinorial Hilbert space. For example, 
we have already discussed this problem in detail in section |^ (there symmetry principles such 
as general coordinate and local Lorentz invariance of the index are central considerations). The 
key observation of the previous subsection in this respect is the role of the coupling to the 
trace of the torsion V^. Needless to say exactly this coupling is absent in the classical Poisson 
bracket formulation of the Killing- Yano technique. Upon inclusion of the torsion trace one must 
reanalyse the quantum anticommutator of Dirac operators. In explicit examples we have found 
that in order to obtain a vanishing anticommutator one must ensure that both p (e, V) and 
p (e. A) are antihermitean on the manifold in question which necessitates the addition of the 
abovementioned torsion trace V^. 

Unfortunately, manifolds which admit Killing- Yano tensors are not nearly as well understood 
as their cousins that play a role in D > 2 extended supersymmetry: Kahler and hyper-Kahler 
manifolds. The most extensive systematic study so far was made in [O, H2|. Though very explicit. 



their analysis only concerns the local geometrical structure of manifolds admitting Killing- Yano 
tensors. Lacking a systematic topological analysis, this information has to be extracted for every 
example being studied. For this reason, we will in this paper focus on the Taub-NUT manifold 
1^ , 53] for all sample calculations. It should however be noted that at least the Kerr metric 



is tractable as well; an extension of those results to include the Kerr-Newman family could be 
used to study a generalisation of our results to include electromagnetic coupling, although there 
are no conceptual problems expected there. The much simpler Schwarzschild metric only admits 
a non-invertible Killing- Yano tensor, which therefore does not lead to a non-singular inverse 
vierbein. 
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4.3 Example: Taub— NUT geometry and its dual 

Both Taub-NUT and dual Taub-NUT possess an 5*0 (3) isometry generated by 

cos (J) 

ki^dn = ki = - sm(j)dm - cos0cot6'9ui + -^—^Om , 

sm 

k2 = cos(j)d[0] - sin cot 6'9r<^i + ^^r^dm , 

sm tl 

as well as the trivial vector 



(4.17) 



ks = di^] . (4.18) 

(there is also a fourth trivial Killing vector k^ = 9[^]). For the computation of the bulk contri- 
butions to the index, as reported in section [4.4. 1| , it is not very important to make this isometry 
manifest and one can use an (almost) diagonal inverse vierbein based on the metrics ( |1.2| ) or 
( |1.3| ). It is, however, crucial to make this symmetry manifest in order to use the t] invariant 
computations for manifolds with boundary of subsection |3.2| . To achieve this one makes use 
of the left-invariant forms S*. These forms satisfy 

Ck^' = {ikd + dik) = (PS^Sj, + 5,PS;) dx" = , (4.19) 

where Ck is the Lie derivative with respect to any one of the three Killing vector fields ki. 
Explicitly, 

E^^ dx^ = T,'^ = cos ip d9 + sin ip sin 9 dcf) , 

= -sinV^d^ + cos^sin0d0, (4.20) 
cos^d^ + dV-, 

and they satisfy dE^ = -(l/2)e^^'=E* A E-^'. The metric of Taub-NUT then becomes 

d.- = dr^ + r(r + 2m) [(E^)^ + (E^)^] + ^ [E^] ^ (4.21) 

whose determinant is given by 



1 2 



2mr(r + 2m) sin6' (4.22) 



The Killing- Yano tensor can also be expressed using the invariant one-forms, 

Y = eJe/g{e)p^ dx^ A dx"" = -2m dr A E^ + — (r + 2m) (r + m) E^ A E^ (4.23) 

m 

and is therefore 5*0(3) invariant as well. The dual Taub-NUT manifold, obtained through the 
procedure sketched in the present section, is given by the metric 

d.^ = '-±^ dr- + '±±^ + i^r] + ^ [S3] ^ . (4.24) 
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r + 2m 



_ r + m , _ r + m 2 

^12 —z— ^ , ^13 ^ , 

r + zm r + 2m 

'^24 - -T— ^ > W34 — — 2. . 

r + Im r + Zm 



-R43 = R21 = dr A - ^""""^ S2 A 

(r + 2m)2 (r + 2m)^ ' 

~ ~ d A S'^ 2rm ^ ^2 

(r + 2m)3 (r + 2m)2 



Table 5: Geometrical data of the Taub-NUT manifold. 



with determinant 



2mr{r + 2m) sin 6' 
^2 



(4.25) 



where S = [r + m)/m. 

The connections and curvatures are all summarised in the two tables ^ and Before we 
discuss the Dirac operators, let us make a few remarks about the displayed results. First of all, 
the curvature of the dual manifold is not (anti)self-dual, in contrast to Taub-NUT (note that 
in our conventions R^'^ = — (l/2)e'"**"i?sn)- Its asymptotic geometry is also different: whereas 
Taub-NUT is asymptotically a flat metric, the dual metric has a volume element independent 
of r at large radius; it tends to an infinite tube. 

Of course, we must also specify the coordinate ranges for the angular coordinates 6, and ip. 
It is easy to see that for both Taub-NUT and dual Taub-NUT one has 9 G [0, tt) and e [0, 27r). 
For the coordinate ijj one finds ip E [0, Air) by ensuring the absence of the NUT singularity at 
r = 0. In detail, one can write both the Taub-NUT and dual Taub-NUT metrics on a constant 
(6'0)-slice in the vicinity of r = as ds^ = dx^ + jx'^dip'^ via the coordinate transformation 
= 8mr. Requiring this be the metric of a two-dimensional plane in polar coordinates we find 
the quoted result for ip. 
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4rm2 
r + 2m 



S3, 



m?r{r + 2m) 2 
r + m 



3 _ yjm?'r{r + 2m) , 4 

6 — 2-/ , 6 

r + m 



r + 2m 



dr . 



0)14 



r + m , 



r + 2m 



m 



(r + 2m)2 
m'^ 



(r + 2m) (r + m)^ 



Wl3 

W34 



r + m 9 



r + 2m 
— 2m^ ^3 
(r + 2m)2 



m"" 



(r + 2m)(r + m)' 



^12 



1 r 



3 r + m 
1 r 



3 r + m 

2 r 

3 r + 2m 



-R21 

^31 

Rai 



m 



(r + 2m)2 
m 

(r + 2m)2 

4m^ 
(r + 2m)^ 



dr AS^-2 
dr A 2^+2 
dr A S^- 



(r + ml 



m 



(r + 2m)3 (r + m)2(r + 2m)^ 
(r + m)3 



m" 



(r + 2m)^ (r + m)^(r + 2m)'^ 
2m^r 



(r + m)(r + 2m)^ 



A s^ 



A s^ 
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m^fSr + 5m) , 2rm^ n o 
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(r + 2m)2(r + m)'^ 

m'^(3r + 5m) 
(r + 2m)^(r + m)^ 



dr A + 



(r + 2m)'^ 

2rm? 
(r + 2m)'^ 



t} a 
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4m(r + m) 
(r + 2m)^ 



dr A + 



2r + 2rm — m 



+ 



TIT 



(r + 2m)^ (r + 2m)^(r + m)^ 



Table 6: Geometrical data of the dual Taub-NUT manifold. 

For the Dirac operators we actually need the inverse of the invariant one-forms. We denote 
them by Sj just as in subsection p.2|. They are the vectors commuting with the Killing vectors, 



[fc„S,] = 0, (4.26) 

and in components one finds 

V MA V /A , sinV^ cos^sinV^ 

cosV'^ cos?/'cos^ (4.27) 



sm U sm U 



S3 — "9^ • 
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In the manifestly 50(3) invariant basis the Dirac operator on dual Taub-NUT with Killing- Yano 
torsion is given by 



rdr + 



1 / 2r + 3m 



2r 



2 \ r + 2m (r + m) 



I ) C ('^,34^ I ^,42v^ \ I ^^,14v ^,1234-^^ + 2 



. (4.28) 



It is antihermitean on the dual Taub-NUT metric and anticommutes with the operator 



„ 1 / 2r + 3m 

rOr + - 



2r 



2 \ r + 2m (r + m) 



+ <J 7 Li + 7 ^2 + -^7 ^3-7 



. (4.29) 



The operator p (e, y) is antihermitean on the dual Taub-NUT manifold but has no definite 
hermiticity on Taub-NUT itself where it can be regarded as the Dirac operator coupled to the 
trace of the torsion 7^V^ = —'^^/[{r + m) V], or in other words a purely imaginary abelian vector 
field. In fact the replacement 



dr dr + 



r{r + m) 



(4.30) 



renders the operators p (e, V) and p (e, A) antihermitean on Taub-NUT and in the case of the 
former one obtains p{e), the torsion-free Taub-NUT Dirac operator. Note that in the above 
A = 2m /{r + 2m) , 5* = (r + m) /m and V = ^/ {r + 2m) /r. 

4.4 Example: Index theorems for Taub-NUT and dual Taub-NUT 
4.4.1 Bulk 



Using the geometrical data of the two manifolds we just presented, all terms in ( |3.32| ) can be 
computed and the indices can be compared. Let us start with the bulk part. For the dual 
Taub-NUT manifold there are contributions from all terms in expression ( |2.19| ) except for the 
F A F term since F = dA vanishes exactly for this example. Putting the boundary at a finite 
radius r = we obtain the (rather unenlightening) result 



1 



24.87r2 



M 



R^n A i?"'" - -F A F - 2 ^D^1C>' 



1 rl {rl + 12mr^ + SGr^m^ + 340mV^ + 753r^m'' + 872rfem^ + 408m'') 



12 



(rj, + 2m)^(r;, + niY 



. (4.31) 
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In the limit Vb ^ oo this yields the number 1/12. If one turns the torsion off, the bulk result for 
dual Taub-NUT is then simply 

1 f 1 rl {9rt + 42r3m + 73r,^m^ + SGr.m^ + 24m^) 

Kmn/\^ = T-^ 7 -7^ TTT ] • l^.dij 



24.87r2 Jj^ 12 {n + 2m)^{n + m) 



The analogous computation for Taub-NUT is well known |45, 46 



' ^ R^.AR- = l''^^''^:'''':\'''^'\ (4.33) 



24.8772 Jj^ 12 (rb + 2m)^ 

For a boundary at infinity, this produces the answer 1/12 (or equivalently a Pontrjagin number 
2). For Taub-NUT it is known that this number, which seems to lead to a fractional index. 



in fact gets modified by the boundary contributions [0, |1^, |1^]. Because of the general results 
discussed in the first part of this section, we expect that a similar correction will arise for the 
dual Taub-NUT manifold as well. 

4.4.2 T] invariant 

The boundaries of Taub-NUT and the dual Taub-NUT are so our general discussion for 
the computation of the 77 invariant in subsection J?2 applies. From ( 4.28 ) we obtain the dual 
Taub-NUT boundary operator 

e=r'^'^> ° (4 34) 

" -P{-S,Z)) • ^ ' 



analysed above in 3.2 with 



3='-!^ (4.35) 



r + 2m m 



as well as 



Killing- Yano torsion, 
Z = { r (4.36) 
— torsion-free. 

, 2m 

(the torsion-free case is in fact the value for Z at which the rj invariant can be obtained from the 
one of Taub-NUT by replacing A SX, but we have presented the general formula ( ^.26| ) just in 
case one wants to study the interpolation between the two cases). We then find that //(O) with 
the boundary at finite radius r = for dual Taub-NUT with Killing- Yano torsion is 

1 1 rb{rl + lOmr^ + 65r^m2 + 168m^r5 + 168rbm^ -|- 48m^) 

2^b(0) = in + 2mnn + my ^^'^^^ 
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and for dual Taub-NUT without torsion 

1 



1 (3rb + 4m) 2 



-VB{A=0)iO) 



(4.38) 



v;.^. 12 (rb + 2m)4 ' 
For Taub-NUT itself, from (^4.29|) the boundary Dirac operator is 

-rr'-.a,o,)- 

Again A = 2m /{r + 2m) where this result as well as (^4.34|) are given in the same chiral ba- 
sis ( p.lO| ) to facilitate comparison of the space of boundary solutions with positive eigenvalues in 
subsection [4.4.4| below. Clearly the analysis of subsection ^]2| applies equally well to the Taub- 
NUT boundary Dirac operator (|4.39|) and importantly the solutions still take the form quoted 
in ( |3.18| ). Therefore, for the Taub-NUT manifold, we reproduce the results of and find 



1 rl^Am + Tf,)' 



(4.40) 



2''^'"' 12 {n + 2mY 

Note that the expressions (|4.4CI|) and ( [4.37|) go to —1/12 when the boundary is moved out to 
infinite radius, r;, oo. On the other hand, the limiting value for ^//(O) for dual Taub-NUT 
without torsion is —3/4. 

Finally we note that the number of harmonic boundary spinors h vanishes for Taub-NUT 
and dual Taub-NUT with or without torsion. 

4.4.3 Chern— Simons correction 

For the dual of Taub-NUT, we have computed the difference between the boundary integral at 
infinity for the product metric and the same integral for the full metric to be (the contribution 
A(\.A from the abelian axial vector field drops out since in this example F = dA = 0) 



24.87r2 



dM 



J A - A J A J -2)C^ -uj A R+luj Au Auj + 2)C 



2 rf,m^{ 
3 



-Art 



23rfm 



AOm'^rl 



18rf,m^ + 6m^ 



(4.41) 



{rb + 2m)*{rb + m)^ 

Since the dual metric is already a product metric at infinity, one expects this correction to vanish 
in the large Vb limit, which is indeed the case. Turning off the torsion this reduces to 



1 



24.87r2 



dM 



UJ 



aJ aJ - u A R + lu Auj Au 



4^2 



m r, 



3 {vb + 2m)^(rfe + m)^ 

(4.42) 



Finally, for Taub-NUT itself, we confirm the result of 
1 



which in our coordinate system is 



24.87r2 



dM 



uj^ A R^ - luj^ aJ Auj^ - UJ A R + ^uj Auj Auj 



2^2 



m r 



3 {n + 2my 



(4.43) 



The above expression also vanishes in the large Vb limit. 
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4.4.4 Equality of indices 

Although it is now easy to verify the equality of the indices on the Taub-NUT and dual Taub- 
NUT manifolds, it would be remiss of us to do so without first checking the validity the assump- 
tion of subsection |4.1| above. Therefore we return to the solutions presented in equation (|3.18| ) 



and check that the eigenspaces generated by the solutions with positive eigenvalues coincide for 
Taub-NUT and dual Taub-NUT. Clearly the eigenvalues (setting Z = 0, since Z = r/{2m) ^ 
corresponds to dual Taub-NUT without torsion) p/2A + A/4 and \/ 4: + Apq{\Sy + {p — qY /2\ 
(5 = 1 for Taub-NUT and S = {r + m)/m for dual Taub-NUT) are always positive. Further- 
more, from the form of the explicit solutions, we see that the eigenspaces generated by the 
corresponding eigenvectors are identical for both manifolds. (Strictly, one must take r;, suffi- 
ciently small for the preceding statement to hold. However, this is sufficient since, as we are 
about to show, neither manifold possesses harmonic boundary spinors at any positive r?, so that 
for both manifolds, the index is independent of > 0.) 

If the solutions corresponding to the negative root A/4 — ^J Apq^XSy + ip — qY /2\ were 
all to have strictly negative eigenvalues clearly we would be done since the projective boundary 
conditions then coincide for the two manifolds. Yet for any non-negative radius r to the boundary, 
the eigenvalue with the negative root can never be positive for either manifold. To see this we 
just need to look at the minimum of the square root \J Apq^XSY + (p ~ which occurs at 
p = 1 = q. But for both manifolds A/4 — S* is negative so long as r is non-negative which 
completes our verification of the assumption. 

We are now justified in comparing the indices of Taub-NUT and dual Taub-NUT. There- 
fore, adding the bulk contribution (|4.31|) , the r] invariant ( [4.37|) and the Chern-Simons correc- 



index^ (e. A) 



tion ( [4.41| ) of the dual Taub-NUT manifold, we find that the r;, dependence drops out completely, 

rl + I2rlm + SGr^m^ + 340r^m^ + 753r^m^ + 872r^m^ + 408r^m^ 

12(r6 + 2mY{n + mY 
rl + I2rlm + 86r^m2+ 308r^m3+ 569r^m^+ 552rj^m^+ 264r^m^-|- ASrb'm'^ 

12{rb + 2mY{n + mY 
32r^m^ + 184r^m^ + 320r^m^ + lAArlm^ - ASnm^ 
I2{rb + 2mY{rb + mY 



0. 



(4.44) 

(the 7] invariant has to be added rather than subtracted because our conventions for the Dirac 
matrices imply that the upper component states have negative chirality). This result matches 
precisely the (well known) result for Taub-NUT obtained by adding ( |4.33| ), ( |4.4CI| ) and ( |4.43| ) (it 



is easy to verify that there is no additional contribution at any value of to the Taub-NUT 
index from the abelian coupling V^) and therefore 

indexf (e. A) = = indexp {e, V) . (4.45) 
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It is interesting to also compute the index for dual Taub-NUT without torsion. In that case, 
adding ( |4.32|) , ( ^4.38|) and ( [4.42|) , we again find a vanishing result 



index ^ (e, A = 0) = (4.46) 

at any radius to the boundary r?, (one might have even argued that since the dual metric can be 
smoothly deformed into that of Taub-NUT, this was to be expected, although in the presence 
of boundary this statement should be carefully re-examined). 

Even though the bulk and boundary terms decouple completely in the limit ^ oo for 
both manifolds, this is not true for general values of indicating the non-trivial nature of our 
result. Moreover, the torsion terms in the various contributions yield a nontrivial answer which 



matches precisely with that of the metric terms to produce the ri, independent result ( |4.45| ). This 
concludes our check of the APS index theorem generalised to manifolds with torsion. 



5. Summary and conclusions 

In this paper we have presented the modifications of the bulk and boundary terms in the index 
theorems for manifolds with boundary when torsion is added. In addition we have shown how 
the index theorem for such manifolds can be related when they possess mutually anticommuting 
Dirac operators. All these results can be explicitly verified for the case of the Taub-NUT and 
dual Taub-NUT manifolds and we find that the index for both these manifolds vanishes. 

Many important technical and physical issues were solved en route to the above results. In 
particular the appearance of the Nieh-Yan tensor at order 0{(3~^) in the bulk index computation 
has previously caused some controversy []T5|, |I6|, ^ 0, However, once one studies well- 
posed index problems by carefully imposing boundary terms for spinors, the index should be 
independent of the inverse regulator (mass)^ (3. Therefore one suspects that keeping all regulator 
dependent terms in the APS analysis in the presence of torsion, all l/f3 terms should exactly 
cancel. The importance of this remark for the axial anomaly in quantum field theory should not 
go unnoticed. 

Interestingly enough, we note that the torsion does not lead to a modification of the index 
for the dual Taub-Nut manifold, so that one may wonder whether on general grounds torsion 
should be expected contribute to the index or not. We are not aware of any general argument 
suggesting that torsion cannot contribute to the index of manifolds with boundary but from a 
physical viewpoint, where each of the bulk, boundary and generalised Chern-Simons corrections 
are separately of considerable interest, the existence of such an argument would not diminish the 
importance of our results which we now summarise. 

The relation found between manifolds with anticommuting Dirac operators involves the 
so-called Killing- Yano duality. In other words, one must must study the motion of spinning 
particles in curved space and in this paper we have extended existing studies of classical spinning 
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particle dynamics to the quantum case in backgrounds with torsion. The precise path integral 
quantisation was given in section |] and the resolution of possible ordering ambiguities for the 
Dirac operator by invoking anti-hermiticity was discussed in section ^. As mentioned earlier, our 
work may also be viewed as a stringent check of the precise path integral quantisation scheme 
utilised in this paper. 

The relationship we found between manifolds with mutually anticommuting Dirac operators 
depended on being able to understand and compare the Hilbert space of spinors on manifolds 
with boundary. This analysis was made possible by the work of APS and we have shown that 
it can be equally well applied to the torsion-full case also. In particular we note that in all our 
analysis it was possible to take any finite radius to the boundary. This allows one to study 
the non-compact limit in which ^ oo which of course has some topical significance in present 
day studies of anti de Sitter metrics. The agreement we found between indices at finite r?, is a 
very strong check of our results. 

The generalisations to include torsion that we presented of Hitchin's rj invariant computation 
for squashed metrics and to Gilkey's Chern-Simons non-product metric boundary correction 
were in principle straightforward but are of course an important step if one is to understand 
index theorems in the torsion-full case. We note that it would be desirable also to consider cases 
where the antisymmetric contortion also has non-vanishing components in the direction normal 
to the boundary. One might expect then that analogous results to ours would also hold. One 
might also like to proceed in such a case by searching for other manifolds satisfying the Killing- 
Yano duality relation. To this end we note that the Kerr-Newman metric of a rotating black 
hole is certainly tractable along the lines presented in this paper, but we reserve this physically 
interesting metric for further study. 

As a final note we also observe that the generalised Chern-Simons form C = u!AR — ^uA 
uj A uj — 2 JC found in section ^ (where, presumably one ought include the order Nieh-Yan 
term along with possible higher 0{h) corrections in /C) might represent an interesting three 
dimensional field theory in its own right. But again we leave such developments to the future. 
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A. Conventions and general relativity with torsion 

We work exclusively with manifolds of Euclidean signature and our Dirac matrices 7^ are her- 
mitean and satisfy {7*", 7*} = 2 5*^*. Flat (tangent space) indices are denoted by the lower 
case Roman alphabet whereas curved indices are members of the Greek alphabet. Products of 
Dirac matrices are denoted as 7''i'"''" = j^^^ ■ ■ ■7^"! where we (anti) symmetrise with unit weight 
and 7^ = -^e'^'^*'^'~^rstu = 7^^^^* = 7^^ We will often employ differential form notation in which 
d = dx^d^ and flat SO {4) indices are usually suppressed and understood to be traced over. For 
a useful review of gravity with torsion see |^ . 



Our convention for the spin connection in the presence of torsion is encapsulated in the 
vierbein postulate, 

= d,ej - Q\,e/ + u;,\e/ , (A.l) 
where the contortion K^^^ is given by 

(~)A pA / I ™ T<r^ ( A 0^ 

In general, this differs from the torsion itself, 

T^ = de'^ + u,'^n^e\ or TV = ^''[H (A-3) 

the precise relation being 

K%u = g'^ + T^.A - T,,^) . (A.4) 

The torsion-free spin connection is expressed in terms of the vierbein by 

We define the field strength of the contortion by 

Ff,ypa = 4 d[^Kyp^] . (A.6) 
Our conventions for the Riemann curvature are 

R{g, T)^yp„ = g^x dyVt^ + Vt^^x^^ i^p - (/i ^ z^) , (A. 7) 

and -R™" = dcu"^" + tu'^p A cj^" = ^R^,y"^"dx^ A dx'^ (the field strength F is similarly normalised 
to be F = dy4 = ^Fp,ydx^ A dx"). Using ( |A.1| ) one verifies that 

-R(^g pppfjCjYi^Cn R(^^^ p,umn ■ (■^•'^) 
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Our sign convention for the Ricci tensor is R^^j = R^'' pa. 

We will mostly be dealing with torsion that is fully anti-symmetric, which we denote by 
writing the symbol A^^^, instead of K'^^^,. In this case, the torsion equals minus the contortion, 
so we will not introduce a separate symbol for it. The axial vector obtained by dualising Ap^y is 
denoted A'^ 

= ^upa ) ^pup = ^\\/9 ^pi/px^'^ ) ^^'^^Apup = ^A'^^Ax , F^i, = 2 d[^A,^-j , (A. 10) 

V 9 

from which one obtains 

Fpupa = -plpiV9^upa]xA^) = IVa^pupaDxA^ , (A.ll) 

using ( [A.6| ) and the Schouten identity. For the contractions of the Riemann tensor, it is useful 
to separate the pieces that only involve a Christoffel connection. When the torsion is fully 
ant i- symmetric, this leads to reasonably compact expressions, 

R{g, A)^^pa = R{g) pupa ~\~ 2 D^p^Ay^pa + 

2A%[^A,],A, (A. 12) 

R{g, A)^a = Rig) pa + DxA\a + A'^^Ax.a , (A. 13) 

R{g,A)=R{g) + A^^'Ax.s- (A.14) 

Dp denotes the covariant derivative without torsion. Expressed in terms of the torsion pseudo 
vector we then have (along with analogous formulae for the Ricci and scalar curvatures) 

R{g,A)f,ypa = R{g)pupa - \DypA'*' ^ eu\pa4> 

(A.15) 

+ T& [gp[pA(7]Ay - gi,\^pAa]Ap - gp[pga]uA Ax\ . 

Finally we note that although the spin connection, in all generality has 24 components which 
decompose under SO {A) as a 4 (the trace), a 16 (terms with mixed symmetry) and a 4 (the 
totally antisymmetric piece), the Dirac operator can only couple via the Dirac matrices to the 
two 4's which one may interpret as the coupling to abelian vector and axial vector fields Vp and 
Ap respectively. Therefore the most general Dirac operator is given by 

p (e, V, A) = ^^e/{dp + \u;{e)pstY' + Vp + -^Ap^') 

where the spin connection with torsion is then 

4 _ 1 



LOprs = Uj{e)prs + gC^frK] " g V 5'(e) ep^paC" rC'' sA" . 

and uj{e)prs and ^'(e) are the torsion-less spin connection and metric determinant computed from 
the vierbein e/. We abbreviate p{e,V,A) to p {e,V) = p{e,V,0), p{e,A) = p{e,0,A) or 
even just p , but our intention should always be clear from the context. 
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